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1. Introduction 

The lingering 3 a discrepancy between the measured [jl]] and calculated [^] values of the 
anomalous magnetic moment of the muon raises the possibility of physics beyond the Stan- 
dard Model. However, a clear understanding of the significance of this discrepancy is complicated 
by the fact that the Standard Model calculation uses additional experimental inputs and models 
to account for the hadronic contributions to a^. Additionally, the hadronic corrections give rise 
to the dominant sources of error in the Standard Model computation. Thus a non-perturbative 
first-principles calculation of the hadronic contributions to could clarify the significance of the 
current 3a effect and potentially reduce the error of the Standard Model result. 

The hadronic contributions to are organized as an expansion in the electromagnetic cou- 
pling a. The leading-order contribution is responsible for the largest single source of error in the 
Standard Model prediction for a^. In this work, we calculate the leading-order hadronic correction 
to using two-flavor lattice QCD. We use the gauge field ensembles of the European Twisted 
Mass Collaboration [Qj. By introducing a new method to overcome the difficulties of previous 
calculations, we achieve a statistical precision of less than 2%. All sources of systematic error 
have been examined and appear to be no larger than the statistical error. As a check, we calculate 
the leading-order correction for the electron and tau leptons and reach similar conclusions. The 
final results of our two-flavor calculation, including an estimate of the systematic errors, will be 
presented in [|J], so here we summarize only the recent advances in our calculation. 

2. Leading-order hadronic contribution 

The expression for the leading-order correction to the anomalous magnetic moment a^ vp of a 
lepton / with mass mi was given by Blum ^ as 

a^ = a 2 dQ 2 - i w{Q 2 /m 2 )U R {Q 2 ) (2.1) 
Jo £r 

where n#(<2 2 ) = n(<2 2 ) —11(0) is the renormalized vacuum polarization n(<2 2 ), Q 2 is the Eu- 
clidean momentum and w(Q 2 /mj) is a known function. The computation of H(Q 2 ) is now a 
standard calculation in lattice QCD. The details of our approach were recently described in 

3. Extrapolation and interpolation 

The lattice calculation of Yl{Q 2 ) is restricted to discrete values of Q 2 and must be interpolated 
and extrapolated to Q 2 = to perform the integration in Eq. [2.1[ To describe the high Q 2 region of 
n(<2 2 ), we use 

o 

n h i g h(e 2 )=c+inG 2 £^(G 2 r. 

«=o 

The constant c is necessary because Tl(Q 2 ) by itself is ultraviolet divergent. The expected ln<2 2 
behavior is explicitly accounted for, but the polynomial sum is added to ensure that Ilhigh provides 
a complete basis in which to expand. The low Q 2 region of Yl{Q 2 ) is dominated by the contri- 
butions from the lightest vector-meson states. Therefore, we take the simple tree-level form used 
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Figure 1: Vacuum polarization function, Yl(Q 2 ). 
This is an example of our calculation of Tl(Q 2 ) and 
the interpolation of the entire range of Q 2 that is cal- 
culated on the lattice. 



Figure 2: Leading-order hadronic correction for the 



muon, a^ vp - This is our calculation of a£ v ', using the 
modified method described in the text, compared to 
the two-flavor contribution to the experimental value. 



hvp 



in effective field theories of vector-mesons [0] and then add a polynomial expansion to provide a 
complete basis of functions. The resulting expression is 

n low (G 2 ) = -| 1^755x4- + I "n(Q 2 T • 
9g Q 2 + m 2 Vi ^ 

The sum over i gives the tree-level contribution for M vector mesons with masses ray; and couplings 
gvi- These masses and couplings have precise meanings of their own and are calculated in the same 
lattice computation. The form for Ili ow is then matched to IThigh to provide a complete description 
that is suitable for numerical integration. An example of this interpolation and extrapolation is 
shown in Fig.|l|. This method gives a fully non-perturbative determination of a} vp that does not rely 
on QCD perturbation theory. 



4. Modified method 

We introduce a class of quantities that have the same physical limit as a^ v but approach that 
limit more smoothly as a function of the pseudo-scalar meson mass nips- For any hadronic quantity 
H, we define 

= a 2 1 dQ 2 ^ 2 w(Q 2 /mj ■H 2 hy jH 2 )U R (Q 2 ) 

where H is understood to be calculated at each value of m PS and // p h ys = H(m PS —> m n ). By 
construction a^ v = as mps approaches m n . Any choice of H leads to a valid definition but 
H = my, the lightest vector-meson mass, leads to a mild mps dependence and results in a well- 
controlled lattice calculation of in the physical limit. The results for are shown in Fig. ||. 
The two-flavor contribution to the experimentally measured value of [||, is also shown in 
the plot and we find good agreement between the linearly extrapolated lattice calculation and the 
measured value. As a check of the method, we perform the same calculation and comparison for 
the electron in Fig. || and for the tau in Fig. |] 
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Figure 3: Leading-order hadronic correction for the Figure 4: Leading-order hadronic correction for the 
electron, a^ vp . This is the same as in Fig. ^ but now tau, a!j vp . This is the same as in Fig. ^ but now for the 
for the electron. tau. 



5. Conclusions 

We have calculated the two-flavor contribution to the leading-order hadronic correction to the 
magnetic moments of the three charged leptons in the Standard Model: the electron, muon and tau. 
Combining an extensive study of all the systematics with a modification to the previously used 
method, we can reliably reproduce the two-flavor contribution to the experimentally measured 
values. The next systematic error that needs to be addressed, the omission of the strange and 
charm quark contributions, will be resolved by our ongoing four-flavor lattice calculation. This 
will open the possibility of a completely first-principles determination of the leading-order hadronic 
contribution to the Standard Model values for the lepton magnetic moments. 
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